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ABSTRACT 


We investigate the power of weak measurements in 
the framework of quantum state discrimination. First, 
we define and analyze the notion of weak consecutive 
measurements. Our main result is a convergence the¬ 
orem whereby we demonstrate when and how a set of 
consecutive weak measurements converges to a strong 
measurement. Second, we show that for a small set of 
consecutive weak measurements, long before their con¬ 
vergence, one can separate close states without causing 
their collapse. We thus demonstrate a tradeoff between 
the success probability and the bias of the original vec¬ 
tor towards collapse. Next we use post-selection within 
the Two-State-Vector Formalism and present the non¬ 
linear expansion of the expectation value of the measure¬ 
ment device’s pointer to distinguish between two prede¬ 
termined close vectors. 


INTRODUCTION 


Weak measurement [1] has already been proven to 
be very helpful in several experimental tasks [2-5], as 
well as in revealing fundamental concepts [6-10]. Tasks 
traditionally believed to be self contradictory by nature 
such as determining a particle’s state between two mea¬ 
surements prove to be perfectly possible with the aid of 
this technique. Within the framework of the Two-State- 
Vector Formalism (TSVF), weak measurements reveal 
several new and sometime puzzling phenomena. For a 
general discussion on weak measurements see [1, 11-14]. 
In this paper we analyze the strength of weak measure¬ 
ments by addressing the question of quantum state dis¬ 
crimination. 


discrimination, i.e. for the task of deciding which vector 
was chosen out of a predetermined set of (possibly close) 
state-vectors (for a general review see [15]). Discriminat¬ 
ing between predetermined non-orthogonal vectors can 
be seen as quantum hypothesis testing [15]. Suppose we 
use projective measurements, then the question is what 
are the best projections to choose so that the error prob¬ 
ability would be as small as possible. This was first dis¬ 
cussed by Helstrom in [16] for the case of two predeter¬ 
mined vectors. The density matrix version of the prob¬ 
lem was developed by Osaki in [17]. A different scheme 
was presented by Ivanovic [18] where the discrimination 
is error free but there is a probability for obtaining a non- 
conclusive result, i.e. a ‘don’t know’ result. A variant of 
this scheme was suggested in [19], [20] where the incon¬ 
clusive outcomes have a fixed rate. Recently, Zilberberg 
et al. [21] have presented a new method, based on par¬ 
tial measurements followed by post selections. The prob¬ 
lem of state discrimination and the problem of cloning 
are deeply connected. A recent paper of Yao et al. [22] 
discusses approximate cloning and probabilistic cloning. 
Weak measurements were also proposed for the task of 
discrimination between two very close states [23]. The 
scheme we suggest here in Ch. II is similar, yet somewhat 
more general, while the scheme we use at Ch. I, based 
on an sequential weak measurements is quite different. 
The advantages and disadvantages of each method will 
be discussed. 

In Ch. I we discuss the orbit of a two dimensional 
state-vector under the set of transformations induced by 
weak measurements. This process can be described as a 
biased Gaussian random walk on the unit circle. The 
probability amplitude that govern the next step (the 
‘coin’ probability amplitude) is changing with each step. 
In this, our walk is similar to the one presented in [24], 
but differs from the quantum random walk presented in 
[25]. We thereby rotate two different vectors in oppo¬ 
site directions (therefore in a non-unitary way), where 
we can use a single (strong) projective measurement to 


There are several known methods for quantum state 
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distinguish between them. We show that using enough 
weak measurements (the number of which is a function of 
the weak coupling), the overall success probability con¬ 
verges to the known optimal result for discrimination by 
projective measurements [16]. 

Next we use a different approach; we reduce the num¬ 
ber of measurements, thus compromising the success 
probability, however, gaining an advantage by avoiding 
the collapse. 

In Ch. II we apply the TSVF of weak measurements. 
By choosing the right Hermitian operator and a proper 
post-selection we can get imaginary weak values. Imag¬ 
inary weak values are best suited for the analysis of the 
coordinate variable of the measurement space. We ap¬ 
ply a non-linear expansion of the weak value and use it 
to compute the first and second moments of such a vari¬ 
able. These moments of the coordinate variable change 
as functions of the initial vector. We then pick two state 
vectors maximizing the difference between the two dis¬ 
tributions of the coordinate variable. 


CHAPTER I: A CLONING PROTOCOL USING 
ITERATIVE WEAK MEASUREMENTS 


In this Chapter we shall perform consecutive weak 
measurements (without post-selection) to show that it 
is possible in principle to differentiate between two non- 
orthogonal vectors. Suppose Alice is sending Bob one 
of two predetermined state vectors of a two dimensional 
system S: 


IV'i) = cosc/|0) + sin^jl) 


| V> 2 ) = sin </>|0) + cos cj) 11) 

where |0) and |1) are the eigenvalues of S z and ^ > cj) > 
j. Let 9 be the angle between the two vectors: 


= cos 9, 

therefore the two vectors have the same angle f with 
respect to the vector -^(|0) + |1)) (see Fig. 1). We 


assume Alice is sending each of the vectors with the same 
probability. It is well known by [16] that the maximal 
success probability is: 


Ps(opt) = ^(1 + yjl - 4AiA 2 |(V’i|V , 2)| 2 ) 


where Alice is sending \ipi) (resp. \ip 2 )) with probability 
Ai (resp. A 2 ). Since we are using Ai = A 2 = 1/2 we can 
write: 


1 + sind 2 

Ps(opt) = ---= cos 2 (j) (1) 


Below we shall show that one can reach the same limit 
using the following weak measurement protocol. By a 
series of weak measurements Bob will be able to ‘rotate’ 
the initial vector towards the direction of |0) or |1). Bob 
will stop the rotations after a predetermined number of 
iterations, by then he can assume with high probability 
that the vector would have crossed |0) or |1) which are 
close to the axes |0) or |1). He will then strongly measure 
the final ‘rotated’ vector in the standard S z basis to get 
the result |0) or |1). If the result is |1) he can conclude 
the initial vector was |'0i), otherwise it was 

The error probability has two factors; the first origi¬ 
nates from the weak ‘rotations’, i.e. the probability that 
the weak ‘rotations’ will take |^>i) (resp. 1 ^ 2 )) to |0) 
(resp. |1}). The second factor originates from the strong 
final measurement, i.e. the probability that having ‘ro¬ 
tated’ |'i/ ; i) (resp. 1 ^ 2 )) in the correct direction towards 
|1) (resp. |0» the strong measurements will produce |0) 
(resp. 11)) results. 

Below we start by describing the process of weak mea¬ 
surement. Then we describe the protocol in details. 

1.1. A Gaussian-type random walk induced by 
weak coupling 

Let S denote our two dimensional system to be mea¬ 
sured. Let S z be the Pauli spin matrix on the system S. 
Let | ip) = a|0) + /3|1) be a state vector in the eigenbasis 
|0) and |1) of S z . 

Let | </} denote the wave function of a quantum mea¬ 
surement device. Then, 
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and T is the coupling time. We will use g = 1 throughout 
this Chapter for simplicity. 

Following the weak coupling the system and the mea¬ 
suring device are entangled: 



1) + /3|1) <g> <j)(x + l)]|x)dx 


( 4 ) 


FIG. 1: The choice of axes, vectors and angles that is 
used throughout the paper. 


where the above functions (j){x ±1) are two normal func¬ 
tions with high variance, overlapping each other. 


I <t>) = \<t>( x )) = [ <t>(x)\x)dx 

J X 


( 2 ) 


where X\x) = x\x) is the position operator of the mea¬ 
suring needle. Suppose |</>(x)| 2 is normally distributed 
around 0 with variance er 2 : 


We can write the entangled (unnormalized) state of the 
measured vector and the measurement device as: 



Qs-ir 


a|0) ® |x) + e 


(*+i r 


/3|1) <8> | x)\dx. 


( 5 ) 


We will now strongly measure the needle. Suppose the 
needle collapses to the vector |xo), then our system is 
now in the state: 


4>{x) = (27rcr 2 ) 1,/4 e x 


The function (j>(x) represents the device’s ‘needle’ distri¬ 
bution amplitude. Let P be the momentum conjugate 
operator of the measuring device, such that [X, P\ = ih. 

We shall start the measuring process with the vector: 


\P) ® |0(x)) 


(a=n—1) , , 

[e a|0) 


(*0 + !) 

e ^ /I|l)] ® |x 0 ). 


( 6 ) 


The eigenvalue Xq could be anywhere around -1 or 1, or 
even further away, especially if a is big enough i.e. when 
the measurement is very weak. Note that the collapse 
of the needle biases the system’s vector. However, if cr 
is very large with respect to the difference between the 
eigenvalues of S z then the bias will be very small and 
the resulting system’s vector will be very similar to the 
original vector. 


in the tensor product space of the two systems. We will 
now couple the two systems by the interaction Hamilto¬ 
nian H int \ 


H = Hint = g{t)S z ® P (3) 


Note that the projective measurement on the outer 
needle’s space induces a unitary evolution on the inner 
space of the composite system. Thus we control the evo¬ 
lution of the state by weak measurements. This resem¬ 
bles an adiabatic evolution of a state vector by strong 
measurements (see [26]). 


where g(t) is the coupling function satisfying: 


g(t)dt = g, 


Consider now the orbit of the initial state vector under 
the series of weak measurements. We couple the parti¬ 
cle to the measuring device and then measure the nee¬ 
dle. Next we couple the biased vector to another mea¬ 
suring device and measure its needle. We repeat this 
process over and over again. The measuring needle is re¬ 
calibrated after each measurement, while the particle’s 
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LogNormal fit, n=2.8 a=0.71 Median collapse time vs. a 2 




FIG. 2: Distribution of the number of measurements 
until the collapse, given a fixed cr. 


state accumulates the successive biases one by one. The 
series of biased vectors describes an asymmetric random 
walk on the circle. 

Note that the random walk is continuous and 
‘weighted’ in the sense that the distribution function for 
the next sampling step is changing as a function of the 
location on the circle. Near the axes |0) and |1) it looks 
like a Gaussian random walk. The following protocol 
goes through such a random walk trying to identify the 
point of start. 

1.2. Distinguishing by consecutive iterations of 
weak measurements 

Bob will perform a series of weak measurements to ro¬ 
tate the initial vector | ijj) towards |0) or |1). We use 
a numerical simulation to investigate the random walk. 
The vectors |0) and |1) are very close to |0) and |1) re¬ 
spectively, and they will define the ‘collapse’. 

First we address the task of quantifying the number of 
weak measurements needed to ‘collapse’ the initial vec¬ 
tor as a function of the standard error a of the needle. 
We simulate the probability distribution of the number 
of weak measurements needed to collapse the initial vec¬ 
tor -1=(|0) + |1}) (see Appendix). We start with cr = 20. 
The probability distribution is best fitted ( R 2 > 0.99) 
by a log-normal distribution with p = 2.8 and a = 0.71, 
see Fig. 2. When examining the form of q m in the Ap¬ 
pendix, and applying the Central Limit Theorem under 
weak dependence, the success of the fit turns obvious: 
q m can be described approximately as an exponent of a 
normally distributed random variable and hence behaves 
like a log-normal random variable. 


FIG. 3: Average number of measurements until the 
collapse as a function of cr. 

Next we look at the median of the above results 
(which distribute log-normally) as a function of cr, 
see Fig. 3. By changing cr we observe a quadratic 
relation between the median and cr. In other words, 
Bob needs 0(a 2 ) steps before he knows the vector had 
‘collapsed’ to |0) or |1) with high probability. In this 
simulation we chose |0) = cos(10°)|0) + sin(10°)|l) ; 
|1) = cos(80°)|0)+sin(80°)|l). The result is independent 
of the initial vector. Rotating |0) and |1) toward the 
axes will only multiply the number of steps needed for 
the collapse by a constant factor. 


Second we address the error probability in determining 
the vector’s original identity. Having started with |^i), 
the error probability Err( |t/>i)) is: 


Err{ |Vb» = P w {\<h) |0))P S (|0) |0»+ 

+P w (\i>i) |i))^ s (|i) ->■ |o». 

The error probability for \ipi) is the probability that the 
series of weak measurement iterations will take to |0) 

and the strong measurement will take |0) to 10), plus the 
probability that the weak process will take \ipi) correctly 
to 11), but the strong measurement will take |1) to |0). 

In the next simulation we compute the probability 
P w {\ipi) —> |1)) and P w (\ip 2 ) —> |0)) for a set of ini¬ 
tial vectors |i/>j) where |0) = cosl°|0) + sinl°|l) and 
|1) = cos89°|0) +sin89°|l). Fig. 4 presents the success 
probability as a function of 9. The simulation was per¬ 
formed 1000 times (see the pseudo-code in the appendix) 
where each time we followed the trajectory of the random 
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Probability of Detection vs. Angular Distance. 



FIG. 4: The success probability P w (—> |1)) as a 
function of 8. The solid curve describes the optimal 
success probability for projective measurements, 
P s (\ipi) —> |1)). Note that the success probabilities are 
with respect to different outcome vectors. 


walk until it crossed the boundaries defined by |0) and 
|1). The success probabilities are higher than the best 
separation value in [16]. However, these simulations dis¬ 
regard the possible error in the strong measurements. As 
we increase the angle between |0) and |1), we reduce the 
error of the strong measurement and the success prob¬ 
abilities of the weak process approach the limit in [16] 
from above. 

The following conclusion is only natural: 

Conclusion Err(\ipi)) = sin 2 a. 

We can demonstrate the conjecture by extending the 
angle between |0) and |1). 

1.3. Hypothesis testing with weak measure¬ 
ments 

So far we have used the weak measurements to itera¬ 
tively produce small biases of the initial vector to shift 
it towards one of the axes. The results we got for the 
pointer of the weak measurement apparatus were so far 
ignored. Suppose now we use a very small number of 
weak measurements. We could use the pointers’ readings 
as samples from the vector’s distribution. Moreover, we 
can average over the few values and use the result as 
a statistic. Notice however that we are sampling from 
a distribution that is changing following each pointer’s 
reading. The advantage of such a protocol lies in the 
fact that the vector has not collapsed; if the standard 
error of the weak measurement is large and the number 
of weak measurement is small then the resulting vector 


Prob. of Detection vs. Angular Distance. 
Data is averaged over 5,10,20 weak meas. 



FIG. 5: Success probability for hypothesis testing with 

low number of weak measurements. The solid curve 
describes the optimal success probability for projective 
measurement (see Eq. 1 above). 

is still in the neighborhood of the original one. We will 
show that the distributions of the averages behaves as a 
function of the initial vector and therefore could be used 
to distinguish between the two. 

In the next simulation (Fig. 5) we weakly measured 
the vectors for 5, 10, and 20 times, using a = 3. This 
simulation was performed 5000 times for several different 
angles 8 (as in Fig. 1). It is expected that the average 
value (of weak measurements) for \ipi) (resp. \ 1 jj 2 ) should 
be below 0 (resp. above 0). The success probability was 
computed by the number of times the average value did 
not cross 0. 

Fig. 6 describes the cumulative distribution function for 
the above average values (denoted by x) for the initial 
angle 8 = 50° and initial vector l'b’ 2 ). It can be seen 
that the graphs are similar to a shifted cumulative distri¬ 
bution function of a normal random variable. Moreover, 
the median is fixed regardless of the number of weak mea¬ 
surements. This median is supposed to coincide with the 
theoretical ‘strong’ one when a becomes lower. 

In practice, if we use such a protocol to distinguish be¬ 
tween two vectors (separated by an angle 8 , as in Fig. 1) 
we should be aware of two types of errors. The first comes 
from the fact that the vector is changing with each weak 
measurement. It could possibly drift following the first 
measurement towards the other vector, thereafter stay¬ 
ing in that neighborhood for long. Thus our readings 
will identify the wrong vector. The second type of error 
comes from hypothesis testing. The sample of the aver¬ 
age could be very close to 0 making the decision tougher. 
The simulation above does not distinguish between the 
two types of errors. It only matches the true vector with 
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CDF of (5,10, 20) weak measurements 



X. In II.2 we show how to pick two initial vectors 
that maximizes the difference between the correspond¬ 
ing distributions of Xf in . Our derivation is based on [27]. 


II. 1. The non-linear approximation of weak val¬ 
ues 

Suppose A is an Hermitian operator on the principle 
system S. Let | 0) denote a state vector for that system. 
Let H = 1 . We will also assume A 2 = 1 , this will make 
it easy to write er lf,AX ! h as a power series. Assume also 
that: 


FIG. 6: Cumulative distribution function for the 
average of weak values for the initial angle 9 = 50° and 
initial vector 1 02 ). 


( A) w 


(i>f in |A|0i n ) _ 

(ipfin | Ip in) 


the average value of weak measurements, to give an over¬ 
all success probability. 


Let | 0) denote the wave function of a quantum mea¬ 
surement device. Then, 


CHAPTER II. STATE DISCRIMINATION WITH 
POST-SELECTION: THE NON-LINEAR 
COMPUTATION 


In this Chapter we will use the TSVF of weak measure¬ 
ments to perform a non-linear analysis of quantum state 
discrimination for a general coupling strength. We will 
show that for the rare events where the post-selection 
is successful we have a high probability to identify the 
vector. The motivation for utilizing post-selection was 
already verified in many precision measurements [2-4], 
hence it is natural to examine it for this application. 

In the pre- and post-selected set of weak measure¬ 
ment (see below) we will look at the distribution of the 
pointer’s coordinate variable Xfi n . The moments of Xfi n 
are functions of the pre-selected state | 0i„), the post- se¬ 
lected state 1'0/m), and the operator A. We will show 
that it is possible to pick two initial vectors |0J n ) such 
that the corresponding distributions of Xf in are easily 
distinguished (possibly by one sample). In particular, for 
one of the initial vectors Xfi n will be distributed around 
0 with standard error a while for the other initial vector 
Xfin will be distributed around a with very low standard 
error (almost 0). This will make it easy to differentiate 
between the two cases. 


|0) = |0(x)) = f (j>(x)\x)dx (7) 


where X\x) = x\x) is the position operator of the measur¬ 
ing needle. We will also assume that |0(a;)| 2 is normally 
distributed around 0 with variance a 2 : 


4>{x) = (2na 2 ) 1 ^ 4 e x 

The function 0(x) represents the device’s ‘needle’ am¬ 
plitude distribution. We will now couple the principle 
system and the measurement system by the interaction 
Hamiltonian Hint'¬ 


ll = H int = g(t)A ® X (8) 

(we used the operator X instead of P since we need imag¬ 
inary weak values [28]). Here g(t) is a coupling function 
satisfying: 


In II. 1 we compute the expansion of the weak 
value (0/m|e _ISj42<: ^|0m) using all terms in A and 


git)dt = g, 
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where T is the coupling time. 

Consider now the average bias of the needle: 

We shall start the measuring process with the vector: 


W) (g> 1 4>(x)) 

/ V\ {&fin\X\$>fin) 

{ )fm ~ ($fin\$fin) 

in the tensor product of the two systems. Then we apply 
the Hamiltonian: 

(X| cos (gX) - i(A) w sm(gX)\ 2 )i n 
(| cos (gX) - i{A) w sin(^i:)| 2 } i „ 

e ~ iAjt ' n \il))®\ cf>(x)). 

Recall (A)„, = ib , let a + = -htL. and a~ = then 

Let 

( X) fin = 

| ®fin(x)) = (^fin\e~ wAjC/h \lpin)\(l>{x)) 

a + (X) in + a~(X cos(2 gX)) in + b(X sin(2 gX)) in 
a+ + a~ (cos(2gX))i n + b(sm(2gX)) in 

be the wave function of the needle following the coupling 
and the post-selection. For an observable M on the nee¬ 
dle’s space ^(a;)), let: 

[29]. Now since the needle is symmetric (normally dis¬ 
tributed) we can write: 

(M). 

( )m m) ’ 

(X cos(2 gX)) in = (sin(2 gX)) in = (X) in = 0. 

Hence 


<X),,„ = b{ X [ (n) 

a+ +a-(cos(2gX)) in 

Since A 2 = 1 we can write {ipfin\e~ l9AX \ipin) as 

We shall now use the parametrization b = cot(^), ( a~ = 
and a+ ~ 2ri^y) therefore: 

°° ( — i n X\ 2n 00 (—‘inX) 2n + 1 

(2nV <*/«!*»> + E f2n + IV 

ivw sm(r])(Xsm(2gX)) in /10 ^ 

WV/in -. , w fr, V\\ ’ ^ ' 

1 - cos(?7)(cos(2 gX)) in 

= (V’/m|^m)[cos(ffX) - sin(gi')] (9) 

Since the needle is Gaussian we can write 


where 


/ t\ (' l l > fin\A\lf)i n ) 

[ >W ('IpfinWin) ' 

(cos(2 gX)) in = e~ 2 ^ 2 

and 


and 
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(Xsm(2gX)) in = 2ga 2 e 2(s<t)2 , 


[30] and therefore 


(X) 


v 

fin 


s\n(g)2ga 2 

e 2{gcr) 2 _ cos ( 7 y) 


(13) 


The above ratio has maximal value at cos(g) = e 2 ^ ga '> 2 
which is: 


a + (X 2 ) in + a (X 2 cos(2 gX)) in + b(X 2 sin(2 gX)) in 
a+ + a~ (cos(2gX)) in + b(sm(2gX)) in 


Since the needle is normally distributed it is easy to 
see that 


(X 2 cos( 2 gX)) in = cr 2 e 2(9<t)2 [1 - 4 g 2 a 2 }, (17) 


[32] and therefore 


{X) 


max 

fin 


2 gcr 2 

y/ e A g °Y - l' 


(14) 


<■ X 2 ) fin = a 2 { 


i+ + a~ e~ 2 ^) 2 [l~Ag 2 a 2 } 
a+ + a~ e~‘ 2 ( grT ') 2 


}• (18) 


[31] If the measurement is weak, i.e. g ■ a <C 1, then 


2 go- 

y/e 4( 9<t ) 2 - 1 

hence 


(X)T°? - <7. (15) 

This will be true for g close to 0. However, if g = 0 then 
(X)} in = 0 for all g ^ 0. 

To sum-up {X)J in is a function of two variables g and 
g. If the two variables are correlated such that cos( 7 y) = 
e - 2 (acr) then for g small enough such that g ■ a <C 1 we 
can get (X)J in to be very close to a. But if we hx g 

(however small) and let g go to 0 we can decrease (X)J in 
to 0. The variable 77 is a function of the weak values of 
A. This means that we can tune g such that for gi the 
expectation value (X)Jl n will be close to a, and for g 2 
the expectation value {X)J 2 m will close to 0. 

To compute the variance of the needle after the 
post-selection note that: 


For small enough ga the value of 1 — 4(? 2 cr 2 will be close 
to 1 and therefore (X 2 )fi n will be close to a 2 . Note that 
(X 2 ) fi n does not depend on the weak value since the 
needle has symmetric distribution. 

II.2. Distinguishing between two non- 

orthogonal vectors 

We can now pick two small angles g\ and g 2 such that 
the difference between (X)^ n and {X)J 2 in is almost cr. 
Since g± and g 2 correspond to two initial vectors, we can 
distinguish between the two vectors by estimating the 
value of (X)J in . In particular, consider 


IV’L) =«il°)-A|i) 

for i= 1 , 2 . Also 


Wfin) = ^=( | 0 ) + | 1 » 

Let A be the Hermitian operator: 



(■ X 2 )fin 


(16) 


Then A 2 = 1. It is easy to see that: 
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Blue: N(|i=0, o=3) ; Green: N(u=3, o=0.3) 


0.5 

Let: _ °- 4 

■o 

Q- 

0.3 

0.2 

r ]i = arcos(e~ 2 ^ ga ' > ). 

o. 


Then we can choose |= ai|0) — Pi\V) such that: 

FIG. 7: Illustration of the normal distribution of the 
pointer’s first moment for each of the two initial vectors. 

ai = 7f (cos T + sm T ) 

The hnal variance of the needle in case |= \ip\ n ) is 



{A)l = i 


Pi 


Oii Pi 


Pi 


i , m . 

—=(cos —— sin — 

V2 2 2’ 


m 


(X 2 ) fin -{X)J 


a 2 — a 2 = 0, 


and therefore 

and therefore the needle is normally distributed around 
a with very low standard error. 

a i + Pi _ cot /^_ ^ The hnal variance of the needle in case | ipi n ) = \tp 2 n ) is 

ai — Pi y 2 


Also, for r] 2 close to 0 we can choose | tpf n ) = a 2 |0) -/? 2 |1) ( X 2 ) fin - {X)J ir ^ « a 2 - 0 = a 2 , 

such that: 


a 2 


I , V2 . . T)2s 
—=(cos — + sm — 
>/2 2 2 ; 


P 2 


1 

71 


, m . m . 

(cosy - Sin-) 


hence 


hence the needle is normally distributed around 0 with 
standard error a. 

We can now easily distinguish between the two alterna¬ 
tive distributions, possibly with a single sample using 
standard hypothesis testing (see for example the illustra¬ 
tion in Fig. 7). 

The success of the protocol depends solely on the post¬ 
selection probability, which is: 


a 2 + P 2 

a 2 — P 2 


,,V2. 

COt(y). 


\{'4 , fin\lpin)\ 


Oii Pi 1 2 



This probability will be low since each of the initial 
Therefore, the difference between (A)^ n and (X)Jl n will vectors is almost orthogonal to the final vector, 
be close to a. 
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To sum-up, having post-selected the final vector, the 
probability to correctly guess the right vector could be 
high. Initially, the probability to post-select is low and 
therefore the complexity of the protocol depends mainly 
on the post-selection. 

In the above example, the two initial states were repre¬ 
sented in the same bases. Including an additional linear 
transformation in the scheme we can generalized it to 
the case of two states chosen from two different mutually 
unbiased bases. This scheme might be suitable for quan¬ 
tum cryptography, and indeed, very recently a method 
based on sequential weak measurements was suggested 
for secure key distribution [33]. 


DISCUSSION 


Weak measurement theory challenges some of the most 
basic principles of quantum theory. In a nutshell, it al¬ 
lows the accumulation of information regarding the state- 
vector without forcing its collapse. As we have shown in 
Ch. I, when performed many times, weak measurements 
are equivalent to a single strong one, thus approaching 
the well-known optimal success probability for discrimi¬ 
nation performed by projective measurement. However, 
when performed only a limited number of times, they do 
not collapse the vector but only rotate it. In such a case 
we can still get some weak information about the state by 
reading (collapsing) the needle of the weak measurement 
apparatus. Moreover, when followed by post-selection, 
weak measurements can reveal underlying properties of 
the initial vector, allowing one to perform quantum state 
discrimination in retrospect as discussed in Ch. II. 

The gradual process of ‘getting information while de¬ 
termining the state’ which was demonstrated in Ch. I 
is strictly connected to the old measurement problem 
[34]. By weakly measuring the initial unknown vector 
we slowly collapse it, creating a continuous tradeoff be¬ 
tween our knowledge and its superposition. This process 
can be thought of as a step-by-step decoherence in which 
the measured system ‘leaks’ through a small hole (the 
weak coupling to the measurement device) into the en¬ 
vironment. As opposed to traditional de-coherence this 
process is rigidly controlled and can be stopped at every 
stage, thus enabling much liberty to the experimenter. 
Therefore weak measurement’s important contribution is 
its flexibility. With weak measurement one controls the 
tradeoff between success rate and collapse rate by choos¬ 
ing the strength of the coupling and the number of weak 
measurement. 


APPENDIX: A PSEUDO-CODE DESCRIBING 
THE WEAK ORBIT 


We wish to find the pointer final position q m (the ex¬ 
pectation of its distribution) when performing m succes¬ 
sive weak measurements on a single particle prepared (for 
example) in the initial state -^(|0) + |1)). Let qo be dis¬ 
tributed according to 1V(1 ,<t 2 ) with probability 1/2 and 
according to N(— l,cr 2 ) with probability 1/2. One can 
take u = m> 1. Now let q\ be distributed according to 
IV(1,<7 2 ) with probability: 

_ exp[(qo + 1) 2 /2<t 2 ]_ 

exp[(q 0 + l) 2 /2cr 2 ] + exp[(q 0 - 1) 2 /2<t 2 ] 

and according to N(— l,er 2 ) with probability: 

_ exp[(q 0 - 1) 2 /2<t 2 ]_ 

exp[(qo + l) 2 /2cr 2 ] + exp[(qo — 1) 2 /2<t 2 ] 

Let <72 be distributed according to N(l,cr 2 ) with proba¬ 
bility: 


_ exp[{(q 0 + l) 2 + (<?i + 1) 2 )/2q- 2 ]_ 

exp[((q 0 + l) 2 + (<7i + 1) 2 )/2<j 2 ] + exp[((q 0 - l) 2 + (qi - 1) 2 )/2<j 2 ] 


and according to N(— l,er 2 ) with probability: 


_ exp{((qo - l) 2 + (gi - 1) 2 )/2<t 2 ]_ 

exp[((q 0 + l) 2 + (<?i + 1) 2 )/2 <j 2 ] + exp[((q 0 - l) 2 + (qi - l) 2 )/2cr 2 ]' 


Then q m is distributed according to 1V(1,<7 2 ) with 
probability: 


_ earKEIlo^ + l) 2 )^ 2 ] _ 

exp[( EEo 1 ^ + l) 2 )/2<x 2 ] + ezpKEEo 1 ^ - l) 2 )/2^ 2 ] 

and according to N(— 1,<t 2 ) with probability: 


_ e^KEITo 1 ^-!) 2 )^ 2 ] _ 

exp[{YZo\<H + l) 2 )/2^ 2 ] + exp[( EEo^ - l) 2 )/2<x 2 ]' 
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